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We present the theory for the effects of superconducting pairing fluctuations on the nuclear 
spin-lattice relaxation rate, 1/Tj, and the NMR Knight shift for layered superconductors in high 
magnetic fields. These results can be used to clarify the origin of the pseudogap in high-T c cuprates, 
which has been attributed to spin fluctuations as well as pairing fluctuations. We present theoretical 
results for s-wave and d-wave pairing fluctuations and show that recent experiments in optimally 
doped YBa2Cu307_<s are described by d-wave pairing fluctuationsEQ In addition, we show that the 
orthorhombic distortion in YBa2Cu307_a accounts for an experimentally observed discrepancy be- 
tween 1/Ti obtained by NQR and NMR at low field. We propose an NMR experiment to distinguish 
a fluctuating s-wave order parameter from a fluctuating strongly anisotropic order parameter, which 
may be applied to the system Nd2-^Ce :E Cu04_5 and possibly other layered superconductors. 



PACS numbers: 74.25.Nf, 74.40. +k, 74.25.Ha 
I. INTRODUCTION 

Fluctuations are enhanced in high-T c cuprate super- 
conductors because of their layered structure and their 
small coherence length.tj In contrast to conventional su- 
perconductors, where the transition is very well described 
by a mean field theory, an extended region of one to 
several Kelvin around the transition is expected to be 
dominated by critical fluctuations in the cuprates. In 
this paper we discuss the effects of gaussian dynamical 
fluctuations above T Cl which are observable over a tem- 
perature range T — T c T c , on the nuclear spin-lattice 
relaxation rate and the NMR Knight shift in high-T c su- 
perconductors. For a comprehensive review on the role 
of NMR-NQR spectroscopy in the study of fluctuation 
effects in high-T c superconductors see Rigamonti et al.u 

Pairing fluctuation effects on the spin-lattice relax- 
ation rate have been investigated in the dirty limit for 
static, long-wavelength fluctuations near T c by Kuboki 
and Fukuyama.tj Heym extended these calculations for 
s-wave pairing fluctuations by including the fluctuation 
corrections to the quasiparticle density of statesfl Ana- 
lytic expressions for the static, long-wavelength fluctu- 
ation corrections to the spin-lattice relaxation rate and 
Knight shift were obtained by Randeria and Varlamov 
for ultra-clean and dirty s-wave superconductors.Q We 
extend their calculations to include finite magnetic fields 
and unconventional pairing for general values of the im- 
purity scattering rate. Our calculations and numerical 
results include dynamical fluctuations and short wave- 
length fluctuations summed over all Landau levels. 

Dynamical quantities such as the fluctuation contri- 
bution to the spin lattice relaxation rate, l/2~i, carry 
valuable information on the type of fluctuations and char- 
acteristic scattering rates and lifetimes. Qualitatively dif- 
ferent behavior for the fluctuation contributions to the 
rate is predicted for different symmetries of the order 



parameter fluctuations u Analysis of the fluctuation cor- 
rections to 1/Ti provides information on the elastic and 
inelastic scattering parameters. The sign of the fluctu- 
ation corrections to 1/Ti is sensitive to pair breaking 
and the symmetry of the pairing fluctuations; thus, non- 
magnetic impurities have almost no pairbreaking effect 
on fluctuations with s-wave symmetry, but have strong 
effects for d-wave pairing. In the case of s-wave pairing a 
large positive fluctuation contribution to 1 /T\ originates 
from the anomalous Maki-Thompson (MT) processBEl 
We show that this process is suppressed in zero field 
almost completely for d-wave pairing if the mean free 
path is shorter than 20 coherence lengths, but cannot 
be neglected near the transition in finite magnetic fields 
or in the ultra-clean limit. Fluctuation corrections to 
the quasiparticle density of states (DOS) dominate the 
anomalous Maki-Thompson processes in the case of d- 
wave pairing symmetry for realistic scattering parameters 
in high-T c cuprates. For a recent review on the role of 
pairing fluctuation corrections to the quasiparticle den- 
sitw-pf states in high-T c superconductors see Varlamov et 

aim 

Recent 63 Cu NQR-NMR experiments on optimally 
doped YBCO by Carretta et al.Q were interpreted in 
terms of a pseudogap originating from superconducting 
fluctuations. Other theories for the pseudogap include 
spin-charge separation, preformed pairs, phase fluctua- 
tions, van Hove scenarios. For a recent-jreview of this 
broad topic and references see Randeria.tLil Chubukov et 
al.13 proposed a magnetic mechanism for the pseudogap 
in which "hot" quasiparticles become gapped by a pre= 
cursor spin-density wave. Recent studies by Auler et a/.Ej 
of 63 Cu and 89 Y NMR in YBCO as a function of doping 
were interpreted as evidence for the vanishing of the pseu- 
dogap for "hot" quasiparticles due to antiferromagnetic 
spin fluctuations exactly at optimal doping, whereas a 
pseudogap for "cold" quasiparticles persisted at optimal 
and overdoped samples. Whether the pseudogap is due 
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to pairing fluctuations, spin-density wave fluctuations, 
or more complicated mechanisms may not be easy to de- 
cide, especially in optimally doped materials. The study 
of fluctuation effects in the presence of strong magnetic 
fields may be key to solving this problem. 

Magnetic fields tend to enhance pairing fluctuations 
near the transition temperature as a result— of Landau 
quantization of the orbital motion of pairs.Eil However, 
because the transition temperature is suppressed by a 
magnetic field, pairing fluctuations are typically reduced 
at constant temperature with increasing field. Applica- 
tion of a magnetic field at constant temperature has very 
different effects on the pairing fluctuation contributions 
to 1 jT\ depending on the pairing symmetry. For s-wave 
pairing the rate is reduced with increasing field, whereas 
in d-wave pairing the suppression of the DOS fluctua- 
tions, which have a negative sign, leads to an enhance- 
ment of l/Ti with field! 

In the next section we describe the theoretical frame- 
work for our analysis of fluctuation effects on NMR in 
high magnetic fields in high-T c superconductors. We de- 
rive the fluctuation propagator for a quasi-2D layered 
superconductor and include quasiparticle scattering by 
non-magnetic impurities in addition to pairbreaking by 
inelastic scattering. We incorporate orbital quantiza- 
tion by the magnetic field on the pairing fluctuations 
as well as the effects of d-wave symmetry. In section 
III we discuss the pairing fluctuation corrections to the 
nuclear spin lattice relaxation (NSLR) rate. To lead- 
ing order in T c /Ef the dominant fluctuation corrections 
are determined by Maki-Thompson processes and cor- 
rections to the quasiparticle density of states. We derive 
expressions for these processes appropriate to 2D fluc- 
tuations in a strong magnetic field and present our re- 
sults for the pairing fluctuation corrections to the NSLR 
rate. The field-dependence of the fluctuations is shown 
to be sensitive to the symmetry of the pairing fluctua- 
tions. In section [V we derive the leading order correc- 



tions to the Pauli spin susceptibility and its contribution 
to the Knight shift. The Knight shift is determined by 
the the long-wavelength spin susceptibility, and in con- 
trast to the NSLR rate the fluctuation corrections to the 
spin susceptibility are not very sensitive to order param- 
eter symmetry or impurity scattering. However, dynam- 
ical fluctuations and orbital quantization lead to signifi- 
cant effects on both the rate and the spin susceptibility 
which are essential for a quantitative understanding of 
the pseudogap behavior in high-T c cuprates. In sections 
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and IV we compare our theoretical results with re- 
cent measurements of the pseudogap in the NSLR, rate 
and the Knight, shift performed by Mitrovic et al.U and 
Bachman et aln on optimally doped YBCO in magnetic 
fields up to 30 Tesla. We show that the pseudogap in 
optimally doped YBCO can be accounted for quantita- 
tively by the theory of 2D pairing fluctuations with d- 
wave symmetry.EJB Finally, we show that incorporating 
orthorhombic anisotropy and the allowed mixing of s- 
wave and d-wave pairing fluctuation channels leads to 



a low-field crossover from predominantly s-wave fluctu- 
ations to predominantly d-wave fluctuations which pro- 
vides a natural explanation for the observed evolution 
from the NQR rate to the low-field (below 2 Tesla) 63 Cu 
NSLR rate on optimally doped YBCO. 



II. PAIR PROPAGATOR FOR 
UNCONVENTIONAL PAIRING 

Fluctuating Cooper pairs are described by a propaga- 
tor, L, which derives from the sum over ladder diagrams 
in the .particle-particle interaction channel as shown in 
Fig. Ffl.li9 Our derivation includes impurity scattering for 
a layered 2D superconductor with an isotropic Fermi sur- 
face and a weak-coupling anisotropic pairing interaction 
g. The generalization to anisotropic 2D and 3D Fermi 
surfaces is straightforward. 



L =® + 




C = x + x 



FIG. 1. Upper row: Diagrammatic representation of 
the Bethe-Salpeter equation for the fluctuation propaga- 
tor, L. The vertex, g, represents the pairing interaction, 
the thick solid lines are quasiparticle Green's functions, 
and the block vertex, C, represents vertex corrections due 
to impurity scattering. Thin double lines symbolize vertex 
factors r](ij)) due to the anisotropy of the pairing inter- 
action. Lower row: Diagrammatic representation of the 
Bethe-Salpeter equation for the impurity vertex correc- 
tions. A thick crossed line stands for the impurity scat- 
tering vertex in the Born approximation. The analytic 
forms for these equations are given in Eqs. (^) and (|^). 

The propagator is a function of the total momentum, 
q, of a pair of interacting quasiparticles, their total exci- 
tation energy, u), and, for anisotropic pairing, their rela- 
tive incoming and outgoing momenta, \ti n ,out- In the fol- 
lowing we use cylindrical coordinates (q, 0, q z ) and write 
q as q = {gcos</>, gsin0, q z }. Pairing fluctuations are 
long-lived only for small lo and q, so that the two par- 
ticles which interact have nearly opposite momenta on 
the Fermi surface, i.e. kj n ~ 2kp.i„ and k out « 2kp j011t . 
We assume a cylindrical Fermi surface of radius hp, in 
which case the momenta on the Fermi surface are given 
by kp = {kp cos'0, kp sini/), k z }. 

The pairing interaction is a function of the momenta 
of the initial and final state of quasiparticles on the Fermi 
surface. We denote the angles between the x-axis (chosen 
as the tetragonal a- axis) and kp and k^ by ip and ip', 
respectively. The pairing interaction, V(tp, t(}'), can be 
expanded in eigenfunctions belonging to the irreducible 
representations of the symmetry group of the crystal. We 



2 



denote the eigenfunction with the largest attractive (pos- 
itive) eigenvalue by rji^j)) and neglect for now the other 
sub-dominant interactions in the expansion of V . Thus, 
we write the pairing interaction in the following form: 



V(ip,ip') = r]{tp) ■ g ■ tj(i//) . 



(1) 



Note that we can neglect the small difference between 
and — q in the pairing interaction, since q ~ l/£o •C 
kp, where £o = hvp /2irkBT c is the coherence length. 

The fluctuation propagator L(w,q) describes dynam- 
ically fluctuating Cooper pairs with a wavelength 2-7r/g 
and a frequency u>. Near T c the typical lifetime of a pair- 
ing fluctuation in the clean limit is 
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where a = 7£(3)/8 w 1.05. We set U = fee = 1 except 
when explicitly noted. 

In the case of strong pair breaking with dephas- 
ing time the prefactor hit/SksT is replaced by r^. 
Spatially small fluctuations decay faster than more ex- 
tended fluctuations. Long-lived fluctuations have typical 

sizes larger than ^ ^J T J } T '-. When the temperature ap- 
proaches T c the importance of long wavelength (q — ► 0), 
quasi-static fluctuations (to — > 0) increases until fluctu- 
ation modes start to overlap in space and time. When 
this happens fluctuation modes interact, which defines 
the critical fluctuation regime. In contrast to conven- 
tional superconductors, where this regime is negligibly 
small, it extends over 1-2 K in layered high- 2^, cuprates 
like YBCO.H Our analysis neglects interactions between 
fluctuation modes and thus excludes the critical regime. 

We include the effects of impurities via the standard 
procedure of averaging over impurity.-»ositions in the 
limit of a long mean free path, I 3> kp 1 Jl3 Impurities lead 
to three different effects: they introduce a finite quasi- 
particlc lifetime via the electron self-energy, they gener- 
ate vertex corrections, V, in the particle-hole channels, 
which have to be included to ensure fundamental con- 
servation laws, and they generate a Cooperon-like mode 
in the particle-particle channel, the impurity vertex C, 
which couples directly to the full pair fluctuation prop- 
agator, L. In the case of d-wave pairing impurities lead 
to pair breaking of the pairing fluctuation modes. We 
will use a short hand notation Q = (uji,q,(j>) for the set 
of arguments related to the pairing channel. The im- 
purity vertex (the cross in Fig. fil) is given in terms 
of the impurity scattering rate in Born approximation, 
a = 1/2t:tNf, 



C(e„, Q) = a + aA (e n , Q)C(e n , Q), 



(3) 



where Ao(e n ,Q) is a momentum-averaged irreducible 
pair susceptibility, defined by the formula (with m = 0) 

/■27T 

A m (e n ,Q) = A m (e ni uji,q,(j)) = N F / dt/j [r](tp)] m 

Jo 



(4) 



Here, £k = e(k) — fx is the quasiparticle dispersion rel- 
ative to the chemical potential. Because q <C kp, we 
approximate £ q _k ~ £k — v f ■ Q- The Matsubara Green's 
functions are given by 
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e» - £k + 5 sign(e„) (7 + ^ 



(5) 



where 1/t^ is the inelastic scattering rate and 1/r is 
the elastic scattering rate. We introduce dimension- 
less scattering parameters, a — h/2irTkBT c and — 
h /27TT0fcsT c . The inelastic scattering rate contributes 
to the quasiparticle scattering, but not to the impurity 
vertex C for the fluctuation propagator. Consequently, 
the lifetime of the pair fluctuation propagator is gov- 
erned by Tij,. Note that both scattering parameters, a 
and cc0, are defined in terms of the renormalized tran- 
sition temperature T c = T c CetpSj^ ) , which is given by 
an Abrikosov-Gorkov formulall3ll3 (see App. |b|), so that 
their values range from zero (for the clean limit) to infin- 
ity (e.g. for the critical pair breaking rate). In high-T c 
cuprates the mean free path.— L is typically of the or- 
der of 3-10 coherence lengths;LD a reasonable estimate is 
(- ~ 5£o> which corresponds to a « 0.2. For the pair 
breaking parameter (or dephasing rate), a^,, one usually 
assumes a much smaller value. For example, compari- 
son between theory and experiment for the c-axis fluc- 
tuation magneto-resistance yields t^T c w iP-|hi YBCO 
and BSCCO, corresponding to w 0.02E3 An esti- 
mate of a$ from inelastic scattering of quasiparticles by 
phonons yields ~ (ffj~) 2 j which at T c in optimally 
doped cuprates is f=a 10~ 2 . However, this weak-coupling 
estimate of inelastic pair breaking may be inappropri- 
ate if the inelastic lifetime is due to strong coupling to 
low-frequency boson modes. Strong coupling or large 
inelastic pair breaking can have a strong effect on the 
pairing fluctuation .eerrections to the nuclear spin lat- 
tice relaxation rateJ13 In weak coupling s-wave theory a 
sign change in the fluctuation corrections to the rate oc- 
curs for « 0.26.E3 A similar sign change occurs in 
strong-coupling theory for a coupling constant A ~ 2. 
Note, however, that a coupling strength of A « 2 is much 
larger than that in conventional strong-coupling super- 
conductors like lead. We consider parameters a<p> 0.26 
and A>2 as unreasonably large for high-T c cuprates. In 
high-Tc materials pair breaking by inelastic scattering is 
probably not strong enough to produce such qualitative 
changes in the behavior of the fluctuation corrections to 
the spin-lattice relaxation rate. Thus, the remaining dis- 
cussion focuses on fluctuations in weak-coupling layered 
superconductors. 

For a single pairing channel in an isotropic quasi 
2D metal the fluctuation propagator factorizes into 



3 



rj(ip)L(Q)i](ip'), where L(Q) obeys the Bethe-Salpeter 
equation: 

L(Q) = g + Tj2gMtn,Q)L(Q) 

n 

+Tj2gMtn, Q)C(e n ,Q)A 1 (e n , Q)L(Q). (6) 

n 

Inserting the Cooperon propagator, C(e n ,Q), from (||) 
into (^|) we can solve for L(Q) in terms of the momen- 
tum integrated pair susceptibilities A m (e n , Q) to obtain, 



HQ) 



1 



where 

B 2 (e n , Q) = A 2 {e n , Q) + A x {e n , Q) 2 C{e n , Q) 



(J) 



A 2 (e n ,Q)+a 



A{ - A A 2 



(e»,Q) 



1 - aA (e n ,Q) 



(8) 



Finally, we must include impurity vertex corrections 
in the particle-particle channel to the external vertices of 
the pair propagator. These corrections are incorporated 
by the replacement rjLrj 1 — ► K with 

K(e n ,e n ,,ip,ip',Q) = rj(e n ,tp,Q) ■ L(Q) ■ fj(e n /,ip',Q), 

(9) 



where 



f)(e n ,ip, Q) = nty) + A! (en, Q)C(e n , Q). (10) 
Combined with Eq. (|^) this gives 

+a(A 1 (e n , Q) - r](Tp)A (e n , Q)) 



rj(e n ,tp,Q) 



1 - aA (e n ,Q) 



(11) 



In the case r](ip) e 1 we recover the standard vertex 
corrections and -pair propagator for an isotropic s-wave 
superconductor ,E2l 



L S (Q) = ■ 



1 - aA Q (e n ,Q)' 
1 



A (e n , Q) 
- aA (e n ,Q) 



(12) 
(13) 



For dynamical quantities such as the spin-lattice re- 
laxation rate it is necessary to analytically continue the 
pair propagator from Matsubara energies to the pal en- 
ergy axis. This is done by Eliashberg's technique,E3 lead- 
ing to the general result 



f T f°° de 



tanh 



e-uj/2 
2T 



tanh 6 j lmB 2 (e,q,(j)) - 2 tanh ^ImB 2c (e) 



+ 



de 



o 2vr v 



( tanh 



e-u)/2 
2T 



, e + ui/2\ 
tanh 2r ) Rc^e, q, < 



(14) 



where 

S 2 (e, g, 0) = r , (15) 



1 - aA (e, q) 
B 2c (e;T) = ^B 2 (^e,q = 0;T = T C ). 



(16) 



Explicit expressions for the functions Aq, A\ and A 2 are 
given for s-wave and d-wave pairing in the Appendix |a|. 

In the long-wavelength limit it is possible to integrate 
Eq. ( |l4| ) analytically and express the pair fluctuation 
propagator for s-wave or d-wave pairing as 



L s {q,io) — Np 1 
L d (q,u)=N F 



1 



1 

(d + £,jq 2 - iujTd ' 



(17) 
(18) 



where the coherence lengths, £ Sj d, static pair susceptibili- 
ties, e S! d and lifetimes, r S)( j are given in terms of Digamma 



functions of the pair breaking parameters (see App. [b|).0 
We generalize the pairing fluctuation theory presented 
above to finite magnetic fields. We assume that the field 
points along the c-axis of the crystal, and introduce the 
following dimcnsionless field, 



b = 



4|e|B 



hvp 



he \27rfcsT c 



(19) 



The main effect of the magnetic field is to quantize the 
orbital motion of the pairs. Through second order in the 
momentum operator, q = — iV — ^fA, quantization of 
the orbital motion is achieved by the replacementslEtEa 



d 2 q 



\b\ ~ 



4 <o to 



(20) 
(21) 
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where k = 0, 1, . . . labels the different Landau levels. 



III. FLUCTUATION CORRECTIONS TO THE 
NUCLEAR SPIN-LATTICE RELAXATION RATE 

The hyperfine interaction between quasiparticles and 
nuclear spins at (fixed) lattice points R„ is given by 

H hf (R u ) = lnle h 2 J d 3 xi(R„) a{x - TL u )S{x) , (22) 

where I is the nuclear spin operator, S(x) is the electron 
spin density operator, and 7 e ,„ are the gyromagnetic ra- 
tios for the electron and nuclear spin, respectively. The 
coupling of the nuclei to the electronic system occurs via 
the hermitian interaction tensor a , which contains the 
contact interaction and dipole-dipole interaction between 
nuclear spin and electronic spin densityH The nuclear 
spin-lattice relaxation rate is well described by second 
order perturbation theory in the hyperfine interaction be- 
tween electrons and nuclei. The transition rate, \j v T™ n , 
from nuclear state | n > to | n' > of a nucleus at lattice 
point Rj, is determined by the matrix elements, "A"™, , 
for the nuclear transition, accompanied by an electronic 
transition from state p to p' , and by the imaginary part 
of the electronic dynamical susceptibility, 



(23) 



-— lim 



dt e i(«+«)t ( [c + a {t}CpP (t) , c+ 7 (0)c k , s (0)]) . 



-'ka 



(cka) creates (annihilates) a conduction electron in 
the Bloch state labeled by k with spin a. We use the 
short hand notation, k = (k, i) for a Bloch state with 
momentum k in band i. The transition rate is given by 



;rpWI 



■ lim 

UJ — »o 



E 

kpotfl k'p'-yS 

Im Xfca,p/3,p' 7 ,fc'<;( a; ) 

UJ 



and all corrections of order (Tc/Ep) 1 arising from pair 
fluctuation modes. We neglect pure weak localization 
corrections and corrections due to the temperature de- 
pendences of the hyperfine coupling matrix elements and 
of the pairing interaction g. Details of the classification 
of diagrams in terms of the small parameter T c /Ep is 
given in Appendix We evaluate the diagrams in the 
Appendix |y. 

The leading order contribution to (TiT) -1 is of or- 
der (T c /Ep)° and defines the Fermi-liquid theory result 
for the normal state NSLR rate, 

{"T-yT)^- = 4n J dk F J dp F N kF N PF rA kFPF | 2 , (26) 

where iVk F is the angle resolved quasiparticle density of 
states on the Fermi surface, and kp defines a point on 
the Fermi surface. The quasiparticle density of states is 
give n by Np = J dkpNk F . The right hand side of Eq. 
( Pq ) is the Korringa constant. til 




(24) 



a) 



FIG. 2. Leading order corrections in T c /Ef to the 
spin-lattice relaxation rate: a) Maki- Thompson, b) and 
c) density of states corrections. V denotes vertex cor- 
rections in the particle-hole channel; V = 1 in our 
model. K denotes the (impurity renormalized) fluctua- 
tion mode in the pairing channel, Eq. (^). 

The fluctuation corrections to X/T\T of order T c /Ep 
are determined in a diagrammatic expansion of the dy- 
namical susceptibility by the Maki-Thompson (MT) dia- 
gram, labeled a) in Fig. 0, and the two density of states 
(DOS) corrections, labeled b) and c) in Fig. ||. The 
Aslamazov-Larkin diagram (not shown) is another order 
smaller in the ratio T c /Ep. The sum of these corrections 
can be written in the following form, 



The matrix elements ^AJJ™ are smooth functions of the 
momenta. Hence, k and p, can be evaluated on the Fermi 
surface. In terms of Bloch wave functions </>/c(x) the hy- 
perfine matrix elements are given by 



v a n'n _ 7"7e fc 2 
A kp — 



h 2 (n\ I{R U ) \ n') 



x J d 3 x 0fe(x) a(x — R„) 4>p(x) 



(25) 



and satisfy "AJJ^ = ( "A™^ n ^ . In what follows we sup- 
press the indices referring to the nuclear transition. 

We perform a systematic expansion of (^TiT) -1 in 
the small parameter T c /Ep (where Ep is the Fermi en- 
ergy) to obtain the leading order term of order (T c /Ep) , 



(TiT)" 1 



E F 



[™ vjqdq [ 2 * 
J 2kN f T c J 



^(S M (q,<t>) + S D (q, </>)). (27) 



where the integrand is-pbtained by analytic continuation 
(following Eliashbcrgca) of the Maki-Thompson and den- 
sity of states correct ions to the dyna mical susceptibility 
obtained from Eqs. ( |D26| ) and ( |D27j ) of Appendix [S], 



J(TiT)" 1 = lim 21m 



Xa/tM + Xdos(u) 



UJ 



(28) 



Thus, we obtain for Sd and S 
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SD(q,4>) 



2 ./o 2tt 



dui i^j — . 
coth— ReL(w,g, 



2tt 



— coth — -ImL (oj, (?, to / — 
2 J 2tt 2T v 'J 2tt 

^ (d w coth ^) ImL(u, q, q 



d l tanh €+ 2 ^ 2 - de tarm 
ft 2 tanh 6 t.^ 2 + ^ 2 tanh 



e-uj/2 

2T 
e-uj/2 



2T 



de 
2^ 



d f tanh — —J- d f tanh 



2T 

e + uj/2 



2T 



-2 / — coth— ReL(^,<j,0) / — 



2tt 



2T 



2tt 



2T 



d f tanh d f tanh 



2T 
e + w/2' 



lm(B N B 2 (e,q, <£)) 

Im (B N B 2 (e,q,<f>)), (29) 



— coth — ImL w.o.ffl / — 



2tt L 2T 

<9 f tanh — h d f tanh 



2tt 



2T 

e + Lu/2 



2T 



f°° duj (■ uj\ , N f°° de 

2 1 ^K coth 2r) Im ^Wo 27 



, e-w/2 
tanh — — tanh 



2T 

e + u/2 



2T 



]m(B 1 (e,q,4>) 2 ) 
Re(Bi(e,g, 0) 2 ) 
Bi(e,g,$| 2 , (30) 



where Bat = 2ttNf, B 2 (e, g, </>) is defined in (15) and 



Bi(e,q, 



Ai(e,q, 4>) 
1 - aA Q (e, q) 



(31) 



The Fermi energy, Ef, is related to measurable proper- 
ties of the 2D Fermi liquid by Ef = fi 2 v F nNFa c , where 
a c is the c-axis dimension of the unit cell. Equation ( j2Sj ) 
originates from corrections to the rate due to pairing fluc- 
tuation corrections to the quasiparticle density of states, 
Fig. H b) and c). The first two terms in Eq. ( p9| ) also 
determine the fluctuation corrections to the Pauli spin 
susceptibility, which we discuss in section IV. Equation 
( ^p| ) represents the Maki-Thompson corrections. The 
first two terms in (^0|) are referred to as the "regular" 
Maki-Thompson contribution, and the last term is the 
"anomalous" Maki-Thompson contribution. The regular 
MT contribution gives a negative correction as does the 
DOS term. The anomalous MT term is positive, but its 
magnitude is very sensitive to pair breaking processes. 
This is the basis for differentiating s-wave and d-wave 
pairing fluctuations using NMR. 

Results for the fluctuation corrections in the quasi- 
static limit are obtained by expanding the integrand for 
small u) (the singularities of the coth-factors are remov- 
able). The long- wavelength limit follows by expanding 
the denominator of the pair propagator to second order 
in q and approximating the remaining terms in the inte- 
grals by their limits for q — > 0. Results for <5(1/T]T) in. 
these limits are discussed by Randeria and Varlamov.Q 
We did not make these approximations; rather we per- 
formed the </>-integral analytically and the integrals over 
e and lo numerically. As we discuss later in this section, 
our approach is important for extending the theory such 
that a quantitative comparison with high-field NMR ex- 
periments can be made. 

In a magnetic field with H\\c, the orbital motion of 
the pairing fluctuations is quantized. Landau level quan- 
tization is achieved by the replacements shown in Eq. 



(p0|). Fluctuation corrections in a magnetic field are of- 
ten treated in the small field limit, where an expansion 
in the magnetic field up to second order is performed. 
At high fields a common approximation is to retain only 
the lowest Landau level. However, one is often in the 
regime between these limits. This is the casepfar the re- 
cent high-field NMR experiments in YBCOMrel To ana- 
lyze this regime we sum the fluctuations over the Landau 
levels numerically. We introduce a cut-off field, b c , to 
regulate the sum over Landau levels, which would oth- 
erwise lead to logarithmically divergent fluctuation cor- 
rections. This divergence is an artifact of the (standard) 
approximation £ q _k ~ 6c — Vf • q made in evaluating 
([|) . Without this approximation convergence is achieved 
on a large momentum scale, ~ hp, or correspondingly 
for n large compared to 1/6. We simulate the conver- 
gence for large n by a cut-off field b c = 20 in our nu- 
merical calculations. Thus, the sum over the Landau 
levels in ( |20| ) extends up to b c /b. Changes in b c lead 
only to overall shifts of the results, indicating small field- 
and temperature-independent "high-energy" corrections. 
These high-energy terms renormalize the leading order 
relaxation rate as discussed below. 



A. Results: Magnetic field dependence 

Calculations of the fluctuation corrections to the 
NSLR rate are shown in Fig. || (for s-wave pairing) 
and Fig. |^ (for d-wave pairing). We normalized the re- 
sults by dividing out the small prefactor, (T C /EA, and 
the normal-state NSLR rate, (TiT)^ 1 ; thus, we plot the 
dimensionless quantity [<5(TiT) _1 / (TiT)^ 1 ] • (E F /T C ). 
Pairing fluctuation corrections in two dimensions contain 
contributions that are constant in temperature and mag- 
netic field. The exact values of these constants are weakly 
dependent on the cut-off in the Landau-level summation 
as mentioned above. These constants, which appear as 
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FIG. 3. Corrections to the spin-lattice relaxation rate 
for T/T c = 95K/92.5K « 1.03 from s-wave pairing fluc- 
tuations as a function of the reduced magnetic field, b. 
The elastic scattering parameter is a = 0.2, and the pair 
breaking parameter, varies as indicated. 
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FIG. 4. Corrections to the spin-lattice relaxation rate 
for T/T c = 95K/92.5K » 1.03 from d-wave pairing fluctu- 
ations as a function of the reduced magnetic field, b. The 
pair breaking parameter is a$ = 0.001, and the elastic 
scattering parameter, a, varies as indicated. 
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FIG. 5. Magnetic field dependence of the Maki— 
Thompson contributions and the DOS contribution 
to the fluctuation corrections to the NSLR rate for 
T/T c = 95K/92.5K ps 1.03, assuming s-wave pairing. For 
comparison, we also show the fluctuation corrections to 
the Pauli spin susceptibility. The curves correspond to 
a = 0.002 and 0.02-0.2 (in steps of 0.02) from top to 
bottom for each set. 




-0.8 -0.4 0.0 0.4 0.8 

b 

FIG. 6. The same as in Fig. |B|, but for d-wave 
pairing; a = 0.002 and 0.02-0.2 (in steps of 0.02) from 
top to bottom for each set. Note, that the anomalous 
Maki-Thompson term dominates for very clean systems, 
a< 0.04, but is negligible for a> 0.1. 



offsets of the curves in all following figures, are irrelevant 
and simply renormalize the normal-state rate, (T\T)~^ . 

Our calculations for the fluctuation corrections to 
1/T\T for s-wave and d-wave pairing symmetry include 
pairbreaking processes from elastic electron-impurity 
scattering and inelastic scattering by emission and ab- 
sorption of phonons. For s-wave symmetry we fixed the 
elastic scattering rate at a = 0.2, and plotted the cor- 
rections for the pair breaking parameter, a^, ranging 
from 0.002 to 0.2. However, for d-wave symmetry non- 
magnetic impurities are already pair breaking, so we fixed 



Q!0 = 0.001 (this value affects the results only in the ultra- 
clean case) and calculated the fluctuation corrections for 
impurity scattering rates ranging from a = 0.002 to 0.2. 
Our results are shown in Figs. || and 1. Note that the 
lowest curve in the d-wave case (Fig. |i) and the highest 
curve in the s-wave case (Fig. 0) correspond to similar 
impurity and inelastic scattering rates, and that s-wave 
and d-wave pairing fluctuations show the opposite field 
evolution in the limit <C a ~ 0.2. Furthermore, the 
s-wave fluctuation corrections to the NSLR rate decrease 
with increasing field even for inelastic rates as large as 
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0(0 ~ 0.1. For very large inelastic rates, > 0.2, the 
maximum in 5(TiT) _1 at b = is suppressed. Such a 
large inelastic pairbreaking parameter appears unlikely 
for the cuprates. More realistic estimates for the elastic 
and inelastic pairbreaking parameters are ~ 0.02 and 
a ~ 0.2OS 

For d-wave pairing the fluctuation correction to the 
NSLR rate changes sign for a w 0.03; the rate de- 
creases with increasing field in the ultra-clean limit and 
increases with increasing field in the limit of weak disor- 
der, a > 0.03. In Fig. ^ we note the rapid drop in the rate 
with increasing field in the ultra-clean limit (a = 0.002) 
compared with the increase in the rate with increasing 
field shown for a = 0.2. It is worth noting that this 
behavior is not obtained in the long-wavelength approx- 
imation employed by other authors.E3 We also note that 
in the clean limit for d-wave pairing the long-wavelength 
approximation is not justified for T ~ Tc > (a + a0).E£l 

In Figs. |^ and || we show the different contributions 
to the relaxation rate for s-wave and d-wave symmetries. 
The label 'DOS' refers to the density of states correc- 
tions in Eq. (p9|). The 'regular Maki' contribution is the 
first two terms in Eq. ( p0| ) and the 'anomalous Maki' 
correction corresponds to the last term in Eq. (|30| ) . The 
full fluctuation correction to ^(TiT)" 1 , shown in Figs. 
H and f|, is the sum of the 'DOS', 'regular Maki', and 
'anomalous Maki' corrections. The 'DOS' term also de- 
termines the fluctuation correction to the tunneling den- 
sity of states at zero bias for an NIS tunnel junction. The 
fluctuation corrections to the spin susceptibility are also 
shown for comparison in Figs. | and § 

For s-wave pairing the regular Maki-Thompson cor- 
rection is (up to a constant) nearly equal to the DOS 
contribution. By contrast, the regular Maki-Thompson 



term is negligible for d-wave pairing. All fluctuation cor- 
rections except the anomalous Maki-Thompson term are 
weakly dependent on the scattering parameters in the 
range of interest. The anomalous Maki-Thompson cor- 
rection is extremely sensitive to pair breaking, as can be 
seen in Figs. || and || Because pair breaking by disorder 
is sensitive to the symmetry of the pairing fluctuations, 
the relative correction to the NSLR rate, shown in Figs. 
^ and [|, shows qualitatively different behavior for s-wave 
and d-wave pairing symmetries. 

In Fig. ^ we show the influence of strong disorder 
on the magnetic field dependence of the NSLR rate for 
an s-wave superconductor. Disorder leads to a reduc- 
tion of the coherence length, and thus to an enhance- 
ment of fluctuations. In the clean limit the typical mag- 
nitude of the fluctuation corrections in 2D contains the 
prefactor, T c /Ep, which is replaced in the dirty limit 
(a = 1/2tttT c > 1) by 1/rEp ~ aT c /E F , which means 
that the fluctuations in dirty s-wave superconductors are 
typically stronger than fluctuations in clean s-wave super- 
conductors with the same T c . By comparison, d-wave su- 
perconductivity is completely suppressed by elastic scat- 
tering for 1/2tttT c q > 0.28, where T c q is the transition 
temperature without impurities. 

Note that the NSLR rate for s-wave pairing decreases 
with increasing field in both the clean and dirty limit 
for realistic pair breaking parameters a0<O.2. The en- 
hancement of fluctuation corrections to the rate reflects 
the reduction in the coherence length by elastic and in- 
elastic scattering. For weak impurity scattering the re- 
duction of the coherence length at T c for s-wave pairing 
becomes 

6(r c ) 2 tt 4 a 

i^= fl -32^ + 3 ) ' (32) 
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FIG. 9. Temperature dependence of fluctuation cor- 
rections to NSLR rate for b = 0.01, a = 0.2 and s-wave 
pairing. 

and for d-wave pairing 

UTcf 



— (a + a), 



(33) 



where a = 7£(3)/8 s» 1.05. Thus, the reduction of the 
coherence length by non-magnetic impurities is stronger 
by a factor of 3 for d-wave pairing compared to s-wave 
pairing at the same T c . This shortening of the coherence 
length is accompanied by a suppression of the transition 
from T c o = T c (l + ^- (a^ + a)) to T c in d-wave symmetry 

compared to T c q = T c (l + ^a^) for s-wave pairing.E3il3 
The slope of b c2 at T c is inversely proportional to the 
square of the coherence length, 



db, 



c2,s(d) 



(IT 



id) 



T=T C 



(d) 



dT 



(34) 



T=T C 



Thus, the reduction in the coherence length leads to a 
significant increase in the slope of b c2 (T) shown in Fig. 
|8|. These results were obtained by numerically solving 
the equation 



, T Ar i f °° de n , e 
In— =Np — 2tanh — 

xlm [B 2 {e,q b ;T)-B 2c {e;T)], 

where q b = VIMW^S 

For weak impurity scattering we obtain 



(35) 



db c2 
dT 



2 
oT c 
1 



« --(1.90 + 0. 
for s-wave pairing, and 



7r / a\ 
Wa {"* + 3 J 

I- 1.83a), 



(36) 
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FIG. 10. Temperature dependence of fluctuation cor- 
rections to NSLR rate for b = 0.01, = 0.001 and d-wave 
pairing. 



db c2 
~~dT 



aT, 
1 

T 



4 ' 32a 

(1.90 + O.81a + 0.81a), 



(a + a) (38) 



(39) 



for d-wave pairing. The negative terms in the brackets 
come from the reduction of the transition temperature 
by pair breaking. 



B. Results: Temperature dependence 

The theory of leading order pairing fluctuations pre- 
dicts characteristic features in the temperature depen- 
dence of the fluctuation corrections to the NSLR rate 
1/T\. Typical results for s-wave and d-wave pairing are 
shown in Figs. |9]-|l6|. For both symmetries there is a pro- 
nounced enhancement of the absolute value of the fluctu- 
ation corrections when the mean-field transition temper- 
ature, T c (b), is approached. However, depending on the 
scattering parameters, a and a^, the corrections may be 
positive or negative near T c (b). 

We first show in Figs. §-|lJ the influence of impu- 
rities on 1/T\ for small and intermediate values of the 
magnetic field, b = 0.01 and b = 0.4. For s-wave symme- 
try we show results for fixed elastic scattering, a = 0.2, 
for a range of pair breaking parameters, a^. In the low- 
field limit, shown in Fig. a crossover from positive 
to negative divergence takes place for w 0.16. The 
divergence is much weaker for strong pair breaking com- 
pared to weak pair breaking in the relatively clean case 
of a < 0.2 discussed in this paper. Note that in the clean 
limit the correction to 1/T\ diverges like y/T c /(T — T c ) 
in zero field,Q .compared to a logarithmic divergence in 
the dirty case.0 In the dirty limit the crossover from a 
positive logarithmic divergence to a .negative logarithmic 
divergence takes place at a^ w 0.26E3 
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The low-field results for d-wave symmetry are shown in 
Fig. [h]. Because inelastic and elastic scattering act simi- 
larly in d-wave superconductors we fixed — 0.001 and 
present results for several values of the elastic scattering 
rate, a. As can be seen in Fig. ^ there is a crossover 
from a positive to a negative divergence for a m 0.04, 
corresponding to a mean free path of about 25 coherence 
lengths. For realistic values of scattering parameters in 
high-T c superconductors, a + » 0.2, a negative diver- 
gence should be observed. 

The effects of a strong field, b = 0.4, are shown in 
Figs. [Il]and[l2| For s-wave fluctuations the pair break- 
ing effect of the magnetic field dominates the effect of 
intrinsic pair breaking, leading to large negative fluctua- 
tion corrections to the NSLR rate as shown in Fig. O. 
For d-wave pairing the effect of a magnetic field is much 
less pronounced. In the clean limit, even at high mag- 
netic fields, the fluctuation contributions to the NSLR 
rate show a positive divergence for d-wave pairing, in 
sharp contrast to s-wave pairing. We discuss this result 
in more detail below. However, for cuprate supercon- 
ductors with d-wave pairing, and a reasonable estimate 
for the scattering rate, a rs 0.2, we obtain a negative 
correction for all field strengths. 

As can be seen by comparison of the NSLR rate for 
b = 0.01 and b — 0.4, there is a strong effect of the 
magnetic field on the temperature dependence in s-wave 
superconductors. The temperature dependences of the 
NSLR rate of superconductors with s- and d-wave pairing 
are compared in Figs. 13 and HJ for different magnetic 



field strengths and parameters a — 0.2, = 0.02_which 
are typical estimates for high-T c superconductors!!-] 

In Fig. [l3| we show, for s-wave pairing, that there 
is a dramatic change in the behavior of the corrections 
to the NSLR rate at field b « 0.2 for T near T c (b). 
Whereas for b<0.2 fluctuations enhance the NSLR rate 
with decreasing temperature, for b> 0.2 fluctuations sup- 



press the NSLR rate with decreasing temperature. Note, 
that to observe this effect one must compare the qualita- 
tive temperature behavior of the NSLR rate for different 
fields rather than changing the magnetic field at constant 
temperature. For d-wave symmetry, shown in Fig. 
this effect is absent. 

To clarify the origin of this behavior we have plotted 
the Maki-Thompson terms, and the DOS term separately 
in Figs. H and ||. As can be seen, all contributions to 
the fluctuations are reduced in magnitude at constant 
temperature with increasing magnetic field. In contrast, 
all terms are enhanced in magnitude with increasing mag- 
netic field for constant T — T c (b), as can be inferred from 
the larger slope of (5(TiT) -1 near T c (b) at lower fields. 

Neither the DOS nor the regular Maki-Thompson cor- 
rection alone are large enough to dominate the anoma- 
lous Maki-Thompson contribution. But together these 
two corrections overcompensate the anomalous MT cor- 
rection for fields above 6>0.2 for s-wave pairing . which 



leads to the qualitative changes shown in Fig. 13, In 
d-wave pairing the regular MT contribution is negligible 
for all magnetic field strengths, as shown in Fig. [16|. This 
is true also for the regular MT contribution in the ultra- 
clean limit, not shown here, and explains why there is 
no change in sign of fluctuation corrections with increas- 
ing magnetic field for d-wave pairing. We also show for 
d-wave pairing in Fig. ^| that the anomalous Maki- 
Thompson term cannot be neglected near T c ; it diverges 
at the mean field transition temperature, T c (b), except 
for zero magnetic field, b. 

Finally, we suggest that the change in sign of the fluc- 
tuation corrections to the NSLR rate for s-wave pairing 
with increasing field should be observable in the electron 
doped compounds like Nd2- a: Ce a: Cu04_5, if they have s- 
wave pairing symmetry. Observation of this effect would 
be a strong confirmation of s-wave pairing in these com- 
pounds. 
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FIG. 13. Temperature dependence of s-wave fluctua- 
tion corrections to NSLR rate for different fields, given as 
the sum of 'anomalous Maki', 'regular Maki' and 'DOS' 
terms. 
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FIG. 14. Temperature dependence of d-wave fluctua- 
tion corrections to NSLR rate for different fields, given as 
the sum of 'anomalous Maki', 'regular Maki' and 'DOS' 
terms. 
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FIG. 15. Temperature dependence of the Maki— 
Thompson contributions and the DOS-contribution to the 
fluctuation corrections in the NSLR rate, assuming s-wave 
pairing. The curves are shown for different fields, ranging 
from 0.01, 0.04-0.4 (in steps of 0.04), from right to left. 
For comparison, we also show the fluctuation corrections 
to the Pauli spin susceptibility. 
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FIG. 16. The same as in Fig. |15[ for d-wave pairing. 
Note, that the regular Maki- Thompson term is negligible 
compared to the other terms at all temperatures. The 
anomalous Maki- Thompson term is negligible at b = for 
a — 0.2 and = 0.02, but contributes considerably at 
higher fields, 6>0.2. 



C. Comparison with experiment 

In order to compare our results with experimental re- 
sults obtained in high-T c cuprate superconductors, we 
discuss first some specific aspects of NMR in these com- 
pounds. In addition to superconducting fluctuations an- 
tiferromagnetic spin fluctuations, are believed to play an 
important role in the cuprates.E£l A spin pseudogap may 
occur at the antiferromagnetic wave vectors q = Qaf, 
which manifests itself in the temperature dependence of 
the NSLR rate of the Cu(2) nuclear spins. The NSLR 
rate is proportional to the slope at zero energy of the 
dynamical susceptibility at the positions of the nuclei, 



i.e. lim^-fO x"(R-v> UJ )/ UJ i an d is especially sensitive to 
changes in the spectral weight of low-energy electronic 
excitations. On the other hand, the Knight-shift ten- 
sor, which probes the static spin susceptibility at q = 0, 
is barely affected by the opening of the antiferromag- 
netic spin pseudogap at Qaf- By contrast the opening 
of a pairing pseudogap at q = affects the quasiparticle 
density of states at the Fermi level, Np, and thus both 
the nuclear spin lattice relaxation rate (~ Np) and the 
Knight shift (~ N F ). 

Recent experiments by Mitrovic et a/.§ and Bachman 
et alu reported the characteristic field scale on which 
the pseudogap behavior is suppressed, H* wlO T in op- 
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timally doped YBCO. Assuming that antiferromagnetic 
correlations lead to a suppression of spectral weight on 
the scale J ~ 100 meV, corresponding to J j hb ~ 1700T, 
this comparatively low magnetic field scale has to be as- 
signed to another origin. Similarly, recent neutron scat- 
tering experiments in fully oxygenated YBCO show that 
the spin fluctuation spectrum near the antiferromagnetic 
wave'vector remains almost unaffected by a field of 11.5 
TeslaEO However, if spin fluctuations are responsible for 
the pairing interaction between quasiparticles, it is possi- 
ble that strong coupling between quasiparticles and spin 
fluctuations may lead to a pseudogap which has char- 
acteristics of both spin fluctuations and pairing fluctua- 
tions. At present a strong-coupling theory of supercon- 
ducting fluctuations and antiferromagnetic spin fluctua- 
tions has not been developed for pseudogap behavior in 
high magnetic fields. 

Our results are based on the theory of weak-coupling 
2D pairing fluctuations. Two-dimensional fluctuation 
theory for YBCO is justified in the presence of magnetic 
fields because of the large vortex liquid region below the 
transition. Phase coherence between planes may be ne- 
glected in the vortex liquid state because of rapid thermal 
motion of the pancake vortices. Thus, it is reasonable 
to neglect the Josephson coupling in the crossover re- 
gion from the normal to vortex liquid state as well. This 
fact, and Landau level quantization in strong magnetic 
fields, implies that fluctuations are predominantly two- 
dimensional. It is possible that for fields smaller than 2 
T a cross-over to three-dimensional behavior might occur 
close to T r . 
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FIG. 17. D-wave calculations for the superconduct- 
ing pairing fluctuation contribution ^TiT) -1 /(TiT)^ 1 of 
63 Cu(2) spin- lattice relaxation rate in optimally doped 
YBCO as a function of magnetic field at temperatures 
ranging from 83 to 102 K ip-i increments of 1 K, and for 
120 K. Circles^ and squaresE3 are NMR and NQR (0 T) 
experiments. The thick curve and the experimental data 
correspond to 95 K. 

We compare our calculations with experiments re- 




cently reported by Mitrovic et alxi on optimally doped 
YBa2Cu30e.95 in a magnetic field H\\c. Our calcu- 
lations, which assume two-dimensional, d-wave pair- 
ing fluctuations describe the experimental NMR data 
remarkably well. The relative fluctuation correction 
5{TiT)~ x / (TiT)^ 1 to 1/TiT for d-wave pairing and a = 
0.2 and several temperatures is shown in Fig. [l7|. We 
define the normal-state rate, (TiT)^ 1 , to include pair- 
ing fluctuation corrections that are constant in temper- 
ature and magnetic field. Thus, to compare with exper- 
iments we subtract these constant shifts from the calcu- 



lated fluctuation corrections as discussed in section III A , 
and define S^T)' 1 = (TiT)" 1 - (TiT)^ 1 . We chose 
the value of the rate at 120K and 30T for this subtrac- 
tion. The experimental results from Mitrovic et aln for 
the fluctuation correction are also shown for the tem- 
perature T = 95K. In order to compare theory and ex- 
periment we subtracted from the experimental data the 
asymptotic normal-state rate, which is well described by 
(TiT)^ 1 ~ T X /(T + T x ), to extract the fluctuation cor- 
rection SiTiT)- 1 . 

The zero-field transition temperature of T c (0) — 92.5 
K determines the absolute temperature scale for the the- 
oretical calculations. We solve numerically Eq. ([35|) for 
the reduction of the mean-field transition temperature as 
a result of Landau quantization. Theoretically the mean- 
field transition temperature is determined by diverging 
pairing fluctuations. To fix the magnetic field scale we 
use the value for the mean field transition temperature at 
8.4 T obtained from our fit to tie fluctuation corrections 
to the Pauli spin susceptibilityj3 discussed in section IV. 



There is one fitting parameter, T c /Ep, which scales the 
magnitude of the fluctuation contributions. As shown in 
Fig. [l7| the agreement between the d-wave fluctuation 
theory and the experimental data from Ref. [l] is excellent. 

We also show in Fig. [l?] data from Y.-Q. Song (black 
squares). cB The data point at H = is the NQR rate. 
The NQR rate is higher than the low field NMR rate in 
the same sample at 3.5 T. A similar drop between the 
NQR rate and the low field NMR rate was obtained by 
Carretta et al. on optimally doped YBCO.B Based on 
the larger NQR rate compared with the NMR rate at 5.9 
T, Carretta et al. concluded that fluctuation .corrections 
to the NSLR rate are predominantly s-waveJa However, 
the field evolution of the NSLR rate from 2 T to 27 T is 
in quantitative agreement with the theory of 2D pairing 
fluctuations with d-wave symmetry, and disagrees quali- 
tatively and quantitatively with the theory of s-wave fluc- 
tuations. The apparent discrepancy between the NQR 
rate and the low-field NMR rate requires explanation. 
We propose an explanation for the low-field evolution 
(0 < ff<2 T) that reconciles Carretta et al.'s sugges- 
tion in terms of s-wave pairing fluctuations with the field 
evolution and our explanation in terms of d-wave pairing 
fluctuations. We show below that sub-dominant s-wave 
fluctuations, induced by the orthorhombic anisotropy of 
YBCO, can account for the low-field evolution. At fields, 
H > 2 T the s-wave fluctuations are suppressed and the 



12 



dominant d-wave fluctuations control the field evolution. 



D. Effect of orthorhombic distortion 



If the crystal symmetry is not perfectly tetragonal, 
then the s-wave and d-wave pairing channels correspond 
to the same irreducible representation. Thus, the pairing 
basisfunction Tj{ip) is of the form 



-4.0 



(40) 



with 01+ 0^ = 1. The results obtained for the fluctuation 
formulae for pure s- and d-wave pairing, Eqs. (|27])-(30), 
are the same with the replacements 



Bi(e,q,(p) = f3 d 
B 2 (e,q,4>)=0, 



+/3 2 A>(e,<7,< 



M{z,q,<l>) a A)(e,g,4>) u ,s 
l-aA (e,q) +Ps l-&A (e,q)' { ) 

2 A 2 {e,q,(j>) + a[A\ - A 2 A ](e,q) 



1 - aA (e, q) 



1 - aA (e,q) 
Ai(e,q,<f>) 



1 - aA (e,q)' 



(42) 
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FIG. 18. Fluctuation corrections to nuclear spin re- 
laxation rate taking into account orthorhombic distortion. 
We assumed an induced asymmetry in the order parame- 
ter described by r)(ip) = /3dT)d(ip) + Psr/s, with 01 varying 
from to 0.2 in steps of 0.025 from bottom to top, and 
d = l-0 s . 



IV. FLUCTUATION CORRECTIONS TO THE 
PAULI SPIN SUSCEPTIBILITY 



Because A\ ~ cos 2<fi the mixed terms in B 2 and B\ which 
enter in the equations ( |29| ) and (^0|) are canceled to a 
large extent by averaging over (f>. So near T c it is a good 
approximation to add the s- and d-wave components of 
the fluctuation corrections with weights 01 and 0^, re- 
spectively. 

In Fig. m we show the result for S^T)- 1 with 01 
ranging from to 0.2. Thus, the low field anomaly in 
the experimental data of Fig. [Tt] can be accounted for 
by a small s-wave component induced by an orthorhom- 
bic distortion, as can be seen comparing with Fig. [lj. 
We estimate 01 ps 0.15 {0\ « 0.85) for optimally doped 
YBCO. Note, that b = 0.8 corresponds to H = 29 T and 
that we account for both the position of the minimum 
in the NMR rate (at w 2 T) and the difference between 
NQR and low field NMR rates with one fitting parameter 

(&)■ 

In BSCCO this effect should be absent if the dom- 
inant pairing channel has Bi g -symmetry (d x 2_ y 2), be- 
cause in this case the lattice distortion does not induce an 
s-wave, but rather a g-wave component with A 2g symme- 
try, which has fluctuation corrections that respond to dis- 
order and field similarly to the d-wave component. How- 
ever, an s-wave component would be induced if the order 
parameter of BSCCO is predominantly B2 9 -symmetry 
(d xy ). 



The Pauli spin susceptibility is obtained from the 
long wavelength limit of the particle-hole susceptibility 
at uj e = 0, 

Xs = A X] X ' **) (°7<5 ' h) Xpy,p8,k a ,kl3(u e = 0), 

(43) 

where h is a unit vector in direction of the applied field 
and fi e = r y e ?i/2. The Pauli spin susceptibility can be ob- 
tained from the spin part of the measured NMR Knight 
shift by subtraction of the orbital and diamagnetic con- 
tributions. Assuming an isotropic hyperfme matrix ele- 
ment, "Ait, and neglecting anisotropic band structure 
and exchange interaction, the spin shift, K spin , is di- 
rectly proportional to the Pauli spin susceptibility, % s . 
The zeroth-order terms in T c /Ep for the particle-hole re- 
sponse function (at lu c — 0) define the Fermi-liquid result 
for the Pauli spin susceptibility, xn- 

The spin susceptibility can be obtained directly from 
the Matsubara Green's functions without analytic con- 
tinuation because it is an equilibrium quantity. Never- 
theless, it is instructive to write down the expression for 
5x in terms of retarded and advanced Green's functions 
defined on the real energy axis. 
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FIG. 19. First corrections in T/Ep to the Pauli spin 
susceptibility. V denotes vertex corrections in the par- 
ticle-hole channel; V — 1 in our model. K, Kl and K2 
denote the (impurity renormalized) fluctuation modes 
in the pairing channel. For the Pauli spin susceptibility 
Kl and K2 either are a singlet Cooperon or a triplet 
impurity Cooperon in a complementary way. 



spin susceptibility are obtained by the procedure dis- 
cussed in App. |Cj, and are summarized by the DOS, 
Maki-Thompson and Aslamazov-Larkin diagrams, shown 
in Fig. [19].c£l Note that in contrast to the large-g e re- 
sponse the contribution d) in Fig. ^ has the same order 
in T c /Ep as the DOS and MT contributions. However, 
it contains only one singlet pair fluctuation mode, the 
other mode in the particle-particle channel is a triplet im- 
purity Cooperon mode. Algebraic expressions for these 
diagrams are given in Appendi x [d|. The sum of the lead- 
ing order corrections in Eq. ( D28[ ) of App. ^ leads to 
the following expression for the relative fluctuation con- 
tribution to the Pauli spin susceptibility: 



5x_ 

XN 



E F 



Vpqdq 
2ttT c 



2tt 



(44) 



The pairing fluctuation corrections to leading order in 
T c /Ep for the static (ui e = 0) long wavelength (q e — > 0) 



where S(q, 4>) sums the contributions from all diagrams 
shown in Fig. [l^ and is given by 



S{q,((>) =tt J — coth—ReL(uj,q,^) J — 



duo uj 



2n 



2T 



2^ 



9i tanh — ^ df tanh 

9* tanh 6 + 9 2 c tanh 



2T 



2T 

e - lo/2 
2T 



ReB 2 {e, q, < 
ImB 2 (e, q, < 



(45) 



Comparing with Eq. (|29|) one realizes that the fluctu- 
ation correction to the Pauli spin susceptibility is given 
exactly by the first two lines of the density-of-states con- 
tribution to the NSLR rate in Eq. (|29|). This result is 
nontrivial not only because the fluctuation corrections to 
the NSLR rate and spin susceptibility are determined by 
different diagrams, but particularly because the NSLR 
rate is a local response defined by an integral over all 
wavelengths, while the spin susceptibility is a global re- 
sponse obtained from the limit q — > 0. The relative cor- 
rections to the spin susceptibility as a function of the 
magnetic field are shown in Figs. |5| and ^ (denoted by 
'susceptibility'). 



A. Results: Magnetic field and temperature 
dependence 



-10.0 - 




-15.0 



-20.0 



FIG. 20. Temperature dependence of fluctuation 
corrections to Pauli spin susceptibility for different 
fields. 



Unlike the NSLR rate the Pauli spin susceptibility is 
not very sensitive to either impurity scattering or order 
parameter symmetry, as can be seen in Figs. ^| and ^ for 
the magnetic field dependence of the susceptibility. Note 
that the small constant offsets have to be subtracted off 
and are included with the leading order terms as dis- 
cussed in section III A. The temperature dependence 
of the fluctuation corrections to \s f° r s-wave pairing 
is shown in Fig. for different magnetic fields. 



The shift in the divergence reflects the field depen- 
dence of T c (b). In Fig. ^l] wc show for comparison 
the magnetic field dependence of the fluctuation correc- 
tions for s-wave and d-wave symmetry. As can be seen 
in this figure, the fluctuation corrections to \s ar e in- 
sensitive to the order parameter symmetry, at least for 
spin-singlet pairing fluctuations. Thus, the mixing of 
s- and d-wave pairing fluctuations due to orthorhombic 
anisotropy, which has a profound effect on the fluctuation 
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FIG. 21. Field dependence of fluctuation correc- 
tions to Pauli spin susceptibility for T/T c = 95K/92.5K 
~ 1.03, and for s- and d-wave pairing. 

corrections to the NSLR rate at low field, has almost no 
effect on the fluctuation corrections to the spin suscepti- 
bility. 



B. Comparison with experiment 



Knight shift measurements in high magnetic fields 
provide valuable information on the fluctuation contribu- 
tions to the Pauli spin susceptibility. The effect of static 
long-wavelength fluctuations on the Pauli susceptibility, 
in zero field have been calculated in three dimensionsEJ 
and two dimensions!] for e = (T-T c ) /T c < 1. The fluctu- 
ation contribution to the spin susceptibility was found to 
scale as Sx/xn ~ m 2D, and 5\/xn ~ const. + ^/e 
in 3D. For the 2D cas e one ob tains (dSx/dT)- 1 ~ T-T c , 
and (dSx/dT)- 1 ~ y/T -T c for the 3D case. Neither of 
these limiting cases is consistent with the recent data of 
Bachman et aln on optimally doped YBCO shown in 
Fig. ~" 
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These NMR measurements of the Pauli spin 
susceptibility do not show singular behavior near the 
transition. This is .typical for a fluctuation-dominated 
crossover transition.0 For this reason it is preferable to 
treat the mean-field transition temperature, T C (H), as a 
fitting parameter. The mean-field transition temperature 
was determined by analyzing the high-precision measure- 
ments of 17 0(2,3) Knight shift in optimally doped YBCO 
at high magnetic fields .□ 

The curvature shown in Fig. |2^ is not reproduced by 
2D static fluctuations in the low-field limit. Three di- 
mensional fluctuations to not account for the behavior 
because they produce curvature in the opposite direction 
compared to the curves in Fig. [2^. We can describe the 
behavior in Fig. ^2] qualitatively and quantitatively by 
taking into account dynamical fluctuations and orbital 
quantization. The magnetic field is in a range where nei- 
ther the low-field approximation nor the lowest-Landau- 
level approximation is applicable. We perform the sum 
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FIG. 22. Calculations for 2D d-wave pairing fluctua- 
tion corrections to the spin susceptibility for B — T - 
14 T (in steps of 2 T). Shown is the inverse of the deriva- 
tive, (dKspin/dT) -1 . Circles show measurements of the 
Kwght shift of 17 0(2,3) in optimally doped YBCO at 8.4 
Tb Open circles denote points used for the fit to K sp i n (T). 

over the Landau levels and over the dynamical modes nu- 
merically. Orbital quantization is the main source of the 
observed curvature at higher fields. Dynamical fluctua- 
tions produce curvature also for zero field, where orbital 
quantization is absent. 

A quantitative comparison of our calculations with 
the experimental data of Bachman et al.u is shown in 
Fig. ^2]. The fit was performed in the region T > 90 K 
directly on the susceptibility data (open circles). Then 
the inverse of the derivative of the experimental data and 
the theoretical curves were calculated; they are extremely 
sensitive to variations at high temperatures where the 
Pauli susceptibility deviates very little from a constant. 
As can be seen in the Fig. |2^, the agreement is excellent 
even up to temperatures of 102 K. The same fit accounts 
for the data in the non-fitted region (full circles) down to 
85 K. The theoretical mean field temperature was deter- 
mined to be about 81 K at 8.4 T. As we discussed in Fig. 
pl| , mixing of an s-wave contribution due to orthorhombic 
anisotropy in YBCO has little influence on the fluctua- 
tion corrections to the Pauli spin susceptibility. 



V. CONCLUSIONS 

We have calculated the pairing fluctuation corrections 
to the nuclear spin-lattice relaxation rate and to the Pauli 
spin susceptibility in 2D s-wave and d-wave high-T c su- 
perconductors in strong magnetic fields. Our calculations 
include dynamical and short wavelength fluctuations. We 
account qualitatively and quantitatively for recent ex- 
periments performed on optimally doped YBCO solely 
in terms of d-wave pairing fluctuations, assuming rea- 
sonable scattering parameters. We find no necessity to 
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invoke the existence of a spin-density fluctuation pseu- 
dogap. We have shown that incorporating orthorhom- 
bic anisotropy and the allowed mixing of s-wave and 
d-wave pairing fluctuation channels leads to a low-field 
crossover from predominantly s-wave fluctuations to pre- 
dominantly d-wave fluctuations which provides a natural 
explanation for the observed evolution from the NQR 
rate to the low-field (below 2 Tesla) 63 Cu NSLR rate 
on optimally doped YBCO. We suggest that a change 
in sign of the fluctuation corrections to the NSLR rate 
near T C (H) with increasing field should be observable in 
the electron doped compounds like Nd2- x Ce x Cu04_5, if 
they have s-wave pairing symmetry Observation of this 
effect would be a strong confirmation of s-wave pairing 
in these compounds. 
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APPENDIX A: IRREDUCIBLE PAIR 
SUSCEPTIBILITIES FOR D-WAVE SYMMETRY 

In this section we summarize expressions for the £- 
integrated Fermi-surface averages of the product Green's 
functions at real energies for the case of d-wave pairing 
in 2D, i.e. rj{\p) — v / 2cos2t/'. The integrals are related 
to Eq. (|4|) by analytic continuation. For d-wave pairing 
they are, 



Me,?) 



2ttYp 



(Al) 



. , s -2ie- J{v F q) 2 - (2i) 2 r- 
Ai{e, q, 0) = A (e, q) — y III ' V = o V2cos2& 



-2i~e + V(uf<?) 2 - (2e) 2 



A (e,q) 



-2ie-y/(v F q)*-(2i)* 



-2ii+y/(v F q)*-(2ey< 



[A 2 1 -A 2 A ](e,q) = 



-A (e,qf 



1 - 



-2ii-y/(v F q)*-(2ey 



-2i~e+^{v F qy-{2ty 



(A2) 
(A3) 

(A4) 



with e — e + inT c (a + a^). 

For isotropic s-wave pairing fluctuations we have triv- 
ially A (e,q) = A 1 (e,q) = A 2 {e,q) given in (jAl|). 



APPENDIX B: COHERENCE LENGTHS AND 
Tc-REDUCTION 

In the long- wavelength, low-frequency limit the pair 
fluctuation propagator for s- and d-wave symmetry 
becomes^ 



L s (q,u) = Np 1 - 

t 



1 



1 



IUJT S 



e d + i z d q 



2„2 



lUJT d 



(Bl) 
(B2) 



This result is obtained by expanding Eq. (|l4|) for small 
q and to and carrying out the e-integral. We define 
cxq = a + a ( p 1 then, 



2T 



_ iTftl + Sulk) + Zlc^'fl 
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J^2 



2T 



16 

2T I 



Td = 



4ttT 

, T t f 1 
In * - 

T r V 2 



e d = In 



T 

~fc. 



1 



2 



2T 

apTc 
2T 



(B3) 
(B4) 

(B5) 
(B6) 
(B7) 

(B8) 



Note that s- and d-wave results differ by more than the 
replacement of by + a. The relative influence of 
a and on the reduction of the coherence length is 
different for d-wave and s-wave symmetry. 

The reduction of T c by impurity scattering is given by 
the Abrikosov-Gorkov formulasEj 



T c0 \2 
T c0 \2 



2 2 
1 2 2 



(s-wave) (B9) 



= (d-wave). (B10) 
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APPENDIX C: CLASSIFICATION OF 
DIAGRAMS 

The essential feature for our classification scheme of 
diagrams is a separation of energy scales. The low en- 
ergy scale set by the temperature {UbT), the quasiparti- 
cle excitation energy (e), the pair excitation energy (w), 
the scattering rates (h/r, %/t$), etc. should be well 
separated from the characteristic high energy scales of 
the metal, e.g. the Fermi energy (Ep). These ener- 
gies define a formal expansion parameter given by the 
ratio of a typical low energy scale and a typical high en- 
ergy scale, for instance, ksT c /Ep. Alternatively one can 
write the formal expansion parameter in terms of the ra- 
tio of a typical atomic length scale (hp 1 , hvp/Ep, etc.) 
and a typical long wavelength scale (£o = hvp /2i:kBT c , 
I = VFT,£<f, = vft$, etc). We perform a system- 
atic expansion in terms of these parameters, and de- 
rive all leading fluctuation corrections in the framework 
of the Green's functions technique. All diagrams pre- 
sented here are understood as containing renormalized el- 
ements. Thus, low-energy fermion Green's functions are 
quasiparticle-Green's functions, vertices are renormalized 
by high-energy quantities. More detailed descriptions of 
this renormalization procedure are given in Refs. 31 33] . 
We assign for simplicity the order of magnitude small to 
the set of expansion parameters (e.g. small = T c /Ep). 
To estimate the order of magnitude of the diagrams we 
replace the Green's functions for the quasiparticles by 
piecewise constant functions, which are equal to 1/ small 
if both the momenta are located in a narrow shell of 
thickness small around the Fermi surface and the ener- 
gies are small, e < small. The corresponding part of 
phase space is called low- energy region. Outside of this 
phase space area, in the high- energy region, we assign to 
the phase space area a measure of one, and the high- 
energy Green's functions are set equal to one. Analo- 
gously, the low-energy range of a pair fluctuation mode 
consists of small pair excitation energies, Tilo < small, 
and small pair momenta, |q| < small. Performing the 
trivial integrations over the step-like Green's functions 
in the asymptotic limit small — > gives the order of the 
diagram. This is done in the following steps: 

1. Estimate the integrand from the number of quasi- 
particle lines in the diagram, uq, which gives a fac- 
tor small~^K 

2. Labeling of the diagram respecting energy and mo- 
mentum conservation. 

3. Estimate the phase space factors: 

a) Restricting all energies to their low-energy region 
gives a factor small^ nE \ where, He is the number 



of independent internal energies. 

b) Restricting the pair momentum to its low-energy 
region gives a factor small D for every quasiparti- 
cle pair in the fluctuation channel, which is not 
otherwise restricted to the low energy region. The 
physical dimension, £), enters explicitly. 

c) Restricting all remaining fermion momenta to 
their low-energy region is the only nontrivial part 
of the estimate. The number of restrictions, nx, 
gives a factor small^ nK \ Note, that the sum of 
two low-energy momenta is not necessarily in the 
low-energy region again. One needs additional geo- 
metrical restrictions to the angles between the mo- 
menta. 

The leading order corrections in 2D to the NSRL rate 
are determined by the diagrams in Fig. |^ for short- 
wavelength external perturbations, hq e ~ pp, while the 
leading order corrections in 2 D to the spin susceptibil- 
ity are determined by the diagrams in Fig. [l^ for long- 
wavelength external perturbations, hq e <C pp. 

In three dimensions these corrections are another or- 
der higher in small, showing the insignificance of fluctua- 
tions in conventional 3D superconductors. In one dimen- 
sion they are of leading order, signaling the breakdown 
of the quasiparticle picture. 



APPENDIX D: CORRECTIONS TO THE 
PARTICLE-HOLE SUSCEPTIBILITY 

We use a short-hand notation, for the combined 
(bosonic) Matsubara-energy and momentum of the the 
pairing fluctuation mode: Q = (w;,q). Similarly, P = 
(e„,p), P' ee (e„/,p'), Q - P = (wj - e„,q - p), 
Y^ P = TJ2 n J2 p etc. We use the usual Feynman rules 
for evaluating diagrams JtB Although we consider spin- 
singlet, s- or d-wave pairing, both spin-singlet and spin- 
triplet fluctuation channels contribute because of triplet 
impurity Cooperons. We neglect the Zeeman coupling of 
the quasiparticle propagators to the magnetic field. This 
allows us to decompose the vertices and fluctuation prop- 
agator in the particle-particle channel into spin-singlet 
and spin-triplet components, 

T aMS {P,P\Q) = 

T S (P, P', Q) al^ lS + r*(P, P', Q) (o%) a/9 (orr») 74 , (Dl) 
K aPlS {P,P\Q) = 

K S (P, P', Q) al fj a v iS + K\P, P', Q) (A) a ^(^) 7i . (D2) 

The Bethe-Salpeter equation for the fluctuation propa- 
gator, 
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KafrsiP, P', Q) = T a0yS (P, P',Q) + |EE T ^ve(P P", Q) ■ G(P")G(Q - P") ■ K e7iyS (P", P', Q) . (D3) 

e,rj P" 



separates into singlet and triplet channels, 

k b *{p, p', Q) = r s '*(p, p',q) + t-J2 rs,i (^ p "i Q) 

P" 

xG(P")G(Q - P") ■ K S ^(P", P', Q) . (D4) 

Corrections to the NSLR rate are described by the 
diagrams in Fig. H The first diagram was investigated 
by Maki and Thompson,B and the last two diagrams 



XDOs{u m ) 



The term xmt corresponds to the Maki-Thompson 
diagram, a) in Fig. |], and the second term, xdos, to 
the two DOS diagrams, b) and c) in Fig. ||. We use the 
relations A_p_p< = (App/)* and Ap^p = (App/)* to 



represent contributions to the NSLR rate from fluctu- 
ation contributions to the quasiparticle density of states. 
Particle-hole vertex corrections, labeled V in Fig. g, can 
be neglected to leading order in T c /Ep above T c because 
they are all proportional to J G(e n , £k) 2 ~ 0. The 
expressions corresponding to the diagrams in Fig. 
with external Matsubara energy uj mi are then (we use 
W = (uj m ,q')) 



simplify the results. 

The fluctuation corrections to the Pauli spin suscep- 
tibility are obtained from the diagrams in Fig. |l9|, 



= Y E (^«A_p, w _p) (<t 7/3 Ap,p_w) G(P - W) G{P) G(Q -P)G(Q-P + W) K a0jS (P, P-W,Q) 

af3~f8 PQq' 

= -%Y \Ap,P-w\ 2 G(P-W)G(P)G(Q-P)G(Q-P + W) (K S {P,P -W,Q) - K\P,P -W,Q)) , 
PQi' 

(D5) 

= Y E (^ 7 Ap_ w ,p) {<T 7a A P>P - W ) (G(P -W) + G(P + W)) G(Pf G(Q - P) K a(i(ja (P, P, Q) 

a/3 7 PQq' 

= 4 Y |Ap-w,p| 2 (G(P~W) + G(P + W))G{P) 2 G(Q- P) (K S {P, P - W,Q) + 3#*(P, P - W,Q)) . 

(D6) 



Xmt(0) = £ X> ■ oSa) (h • t 7/3 ) G(P) 2 G(Q - P) 2 K a ^ s (P, P, Q) 

a(3~/S PQ 

^-2Y J G(P) 2 G{Q-Pf(K s {pP,Q)-K\P,P,Q)) , (D7) 
PQ 

Xdos(0) = 2 Y Y ( h ' "m) ( h ' °7°) G ( p ) 3 G ^ - p ) K ^p a (P; P, Q) 

afdy PQ 

= AYG{PfG{Q-P)(K s {P,P,Q)+m t (P-P;Q)) , (D8) 
PQ 

Xal(0) - Y E ( h ' ( h ■ a vP) G(Pf G{P') 2 G(Q - P) G(Q - P')K aPlS (P P , Q) K Slv( (P', P, Q) 

a/37<5r)C PP'Q 

= ^Y G(P) 2 G(P') 2 G(Q - P) G(Q - P') (K S (P, P', Q)K t (P', P, Q) + K l (P, P', Q)K S (P' 1 P, Q)) . (D9) 

PP'Q 
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The first term, xmt (0) , corresponds to diagram a) in 
Fig. |l^ (Maki-Thompson) , the second term, xdos(0), to 
diagrams b) and c) (DOS) and the last term, Xal(0), to 
the Aslamazov-Larkin diagram, d) in Fig. [19. Particle- 
hole vertex corrections, labeled V in Fig. |19|, can be 
neglected for similar reasons as in the case of the NSLR 
rate. 

To evaluate momentum integrals we split the p-sum 
into a £ p -integral and a Fermi surface average (we use the 
notation (. . .) p := J dp F n{p F ) ■ ■ • where n(p F ) is the an- 
gle resolved (normalized) density of states at the Fermi 
surface). Thus, 



-< 



(D10) 



The lower limit of the integrals is extended from — // to 
— oo. This approximation induces corrections of order 
T c /Ep which vary on a temperature (and field) scale 
large compared to T c and can be incorporated in the 
asym ptotic normal state behavior as discussed in section 
III A . We use the abbreviations 
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G^PQ) 


= N FJ 


f d£ p G(P) 3 G(Q 
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(D15) 



G 4 ,(P, Q)=N F J d£ p G(P) 2 G(Q - Pf. (D16) 

These expressions can be evaluated by complex integra- 
tion. After ^-integration (as on the left hand sides of the 
above equations) the momenta are confined to the Fermi 
surface. The formal identity 



G[E(Cn)](l) 



1 



nil 



G[E(€n)](P) 

(D17) 



where S(e„) is the self energy for the Green's function, 
implies 



G 



2+ro [£(e n ),£(^-e n )](P,Q) 



to! \<5£(e„) 



G 2 [£(e n ),£(^-e„)](P,Q). (D18) 



Because the functional dependence of G 2 on the self en- 
ergies contains only the combination S(e„) — E(cj; — e„), 



we obtain the relations G^, = 



-2CU and G 3 , = -G 3 . 



In the weak-coupling theory for pair fluctua- 
tions we have K s \p,P',Q) = K(P,P',Q) = 
f)(P,Q)L(Q)fj(P',Q), and we can replace X*(P,P',Q) 
and i4r*(P',P, Q) in the Aslamazov-Larkin diagram by 
G(e„, Q)[8(e n > - e„) - 8{uoi - e n > - e n )]/2. The quanti- 
ties f], L and C are defined in section 0. We neglect 
diagrams containing only impurity interactions (and no 
pairing interaction), which describe pure weak localiza- 
tion effects. Furthermore, we assume that the hyperfme 
matrix elements are isotropic on the Fermi surface. Thus, 
we obtain for the NSLR rate, 



XmtM = -2 • |A | 2 • ]T (fj(P, Q)G 2 (P, Q)) p • (fj(P - W, Q)G 2 (P - W, Q)) p _ q , ■ L{Q), (D19) 

e n ,Q 

XdosM= 4-|A| 2 -^ ( ? 7(^Q) 2 G 3 (P,Q)) P (G 1 (P-VF)) p _ q , -L(Q), (D20) 

and for the Pauli susceptibility, 

Xmt+dos(0) = 8 • ]T ^ QfG^P, Q)) P ■ L(Q), (D21) 



XAL 



(°) = 8 ' E ((V(P,Q)G 3 (P,Q)) P ) 2 -C(e n ,Q)-L(Q). (D22) 



These results can be written more compactly using the identities, 
mQ) 2 jJ^G 2 (P,Q)) = -J—( v (p)f}(P,Q)G 2 (P,Q)) , (D23) 



mQ) 2 -^yG 2 (P,Q)^ + 2(n(P,Q)-^-^G 2 (P,Q)\ C(e n , Q) = *-L(p)fj(P, Q)G 2 (P, Q)\ . (D24) 

(D25) 
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Defining Gi(e„) = {^(P)) , B 1 {e n ,Q) = {f,(PQ)G 2 {P,Q)) and B 2 (e n , Q) = (r)(p)rj(P, Q)Q 2 (P, Q)j we obtain 
for the NSLR rate, 



XmtM = -2- |A| 2 • B 1 (e n ,Q)B 1 (e n -w m ,Q)-L(Q), 



XDOs(u m )= 4-|A| 2 - ^ Gi(e 



<5£(e„ 



■B 2 (en,Q)-i(Q). 



(D26) 
(D27) 



and for the Pauli susceptibility, 



Xmt+dos+al(0) — 4 • } 



7 B 2 (e n ,Q)-L(Q). 



(D28) 
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